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The method of the quasiclassical Green's function is used to determine the equilibrium properties of one- 
dimensional (ID) interacting Fermi systems, in particular, the bulk and the local (near a hard wall) density of 
states. While this is a novel approach to ID systems, our findings do agree with standard results for Luttinger 
liquids obtained with the bosonization method. Analogies to the so-called P{E) theory of tunneling through 
ultrasmall junctions are pointed out and are exploited. Further applications of the Green's function method 
for ID systems are discussed. 



1 Introduction 

Quasiclassical methods, developed several years ago and used successfully to describe non-equilibrium 
states in superconductors and superfluids (see, e.g., HI El [51 SJ), have recently been extended to meso- and 
nanoscopic systems. This became possible after the formulation of boundary conditions, including hybrid 
structures and interface roughness ||5l|6l|7l. In recent developments we focussed on spin-effects and spatial 
confinement, e.g., concerning the spin-Hall effect in a 2D electron gas |8|, and spin relaxation in narrow 
wires in the presence of spin-orbit coupling [91. Within the quasiclassical approach, it is also possible to 
study the influence of disorder and Coulomb interaction on the same footing I lOlfTTl, which is manageable 
since the method works on an intermediate level: microscopic details of the wave-functions, on the scale of 
the interatomic distance, are integrated out, leaving equations of motion which can be solved with sufficient 
accuracy. 

As an illustrative example, we apply in this paper the quasiclassical method to determine the equilibrium 
Green's function, and hence the density of states (DoS), of the (spinless) ID model known as Luttinger 
liquid lfT2l [T3l [T4l . This model contains two fermionic branches, linearized near the Fermi points, with 
vf and A^o = (tt^f)"^ the bare Fermi velocity and DoS. Only interaction processes with small (^ kp) 
momentum transfer are considered: standard parameters are 174 and 52, for scattering processes involving 
only one or both branches, respectively. The dimensionless quantities 74 = gi/2nvF and 72 = g2/'2,nvF 
are useful. The spectrum of charge fluctuations is linear, LLi{q) = v\q\; the renormalized velocity, v, and the 
parameter K, given by 



^ = «^ [(1+ 74)2 _^|] ^ ^ 



1+74-72 



1+74+72 



1/2 



(1) 



are often called Luttinger liquid parameters. 

In the following Sect. 2 we discuss the solution of the equation of motion for the quasiclassical Green's 
function in the presence of a fluctuating potential, the latter representing the fermion-fermion interaction. 
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The Keldysh technique is employed throughout. Considering the average - with respect to the fluctuating 
field - of the Green's function and an analogy to the P{E) theory, the field fluctuations are related to 
an effective impedance of the Luttinger liquid, and finally to the DoS (Sect. 3). In Sect. 4 we consider 
an insulating boundary, i.e. a reflecting wall, and determine the DoS boundary exponent. A concluding 
discussion is given in Sect. 5. 



2 Quasiclassical equation of motion and its solution 

In essence, the Keldysh technique ifTSl [161 differs from the zero-temperature and the Matsubara approach 
by employing time-ordering along a contour which runs from —oo to +oo and back to —oo. Rewriting 
the theory in terms of standard (— oo... + oo) integrations, a 2x2 matrix structure results; e.g., the Green's 
function can be cast into the form 

Accordingly, considering some general 2-particle interaction Vq and decoupling this interaction using the 
Hubbard-Stratonovich transformation within the path-integral formalism, two auxiliary fields are neces- 
sary. As a result, the Green's function can be expressed as an average of the free-particle Green's function 
in the presence of fluctuating potentials as follows (see, e.g., Ref. LlOil ): 

G - (Go[$] e^-i"(l+'5o*))o_^ = {Gom)i (3) 

Here Go is flie free-particle Green's function, and Gq ^[$] = Gq^ + l>, where ^ — (j)iaQ + (j)2d'i; (Tq and 
(Ti denote the unit and the first Pauli matrix, and 0i and 02 correspond to the two potentials mentioned 
above. Furthermore, (...)o,$ denotes an average which is Gaussian by construction, and defined such that 

where for clarity space (x) and time (t) arguments are included. The average (...)$, defined through 
Eqs. (O and (|4|i, can be non-Gaussian. (The spin will be suppressed througout this paper) Note that a 
"real" electrical potential can be included by the replacement ipi ^ (f>i + e^cxt in Go[*I'] of Eq. (O; the 
particle's charge is — e. In the next step, we consider the quasiclassical approximation, i.e. we consider the 
difference Gq^ [$]Go [$] — Go [^]Gq^[^] = 0, keeping only the leading terms in the gradients with respect 
to the (spatial) center-of-mass coordinate (which we again denote by x). The equation is then integrated 
with respect to the magnitude of the momentum. The result is 



gtt'ipx) ^i[^,g] (5) 



where p is the direction of the center-of-mass momentum, and 

[l>,g] = ^{xt)gtt'ipx) - gu'{px)^{xt') . (6) 
The quantity gu' {px) is the so-called quasiclassical Green's function, 

gtt'ipx) = ^ J d£_Gu'{p,x), S, = p"^ /2m - p,; (7) 

a major advantage compared to the full Green's function is the normalization, gg — 1. It should be noted 
that in the quasiclassical approximation, only long-wavelength contibutions of the fluctuating fields are 
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taken into account - which is adequate for the problem addressed below. In the ID case, V^; dx, and 
p = ±. For the homogeneous equilibrium case, the solution is 

«(.) = { ; ^!<'> } 

where the Fourier transformed (t ~ t' e) quantity is given, and F{e) — tanh(/3e/2), f3 — l/fc^T. A 
formal solution of Eq. (|5]l is found with the ansatz 

g„Kpa;)=e'^("*)5o,t-t'e-''^(^*') (9) 

provided 

(dt + VFpdxMxt) = DxMxt) = $(xt) . (10) 

Apparently, ip has the same matrix structure as $, namely (p — (pio-Q + tf2^i- For example, we obtain 
explicitly from (|9]l: 

g^{tt') ^ 6{t-t') -2iF{tt')co&ip2Smip'2Sf^-'^'^^ (11) 



g^(tt') = 2F{tt') cos (f2 cos ip'^ e^^'^^^'^'i) (12) 

where Lpj = fj{xt) and tp'^ = (pj{xt') for brevity. We will discuss below that, for the Luttinger model, 
it is sufficient to assume that the phases are Gaussian distributed. It is then straightforward to determine 
(g)^; for example, we obtaii{l| 

- (5(<-<') + ^(«')e"^''^"'Msinli(Ji(tt')) -sinli(J2(tt'))] (13) 

= -{g^Um (14) 

where 

Mtt')^ {{v,-^\f)^/2 = J^"(0) - J^(tt') , (15) 

Ji{tt')= {{^i-p'^){p2+v'2))'f =J'^{tt')-J'^{tt'), (16) 

J2{tt')= ((¥^1 -</p'i)((^2-<^^))* =-j'^{tt')-J^{tt'). (17) 

Note that Jq and J2 are even, and F and Ji are odd under time reversal, t — t' ^ t' — t. Then we define 
J = — Jo + Ji , such that 

Jit) = / ^ - J^iu;)] [B{u;) + 1] (e^ * - l) (18) 

where B{uj) coth(/3tj/2), and we used J^{lu) = [J^(a;) - J^(w)]B(w). From the relation J^(tt') = 
J^{t't) we obtain J^{u}) — J^(— w), implying that J^{uj) — J^{i^) is odd in frequency. Combining the 
above relations, the (normalized) DoS is given by 



N{e) = Re5^(e) = ^ + J dte'"* {F{tt')[P{tt') - P{t't)]}t'=o 

= 1 + ]^ j dujF{e-uj)[P{uj)^P{-Lj)] (19) 



The quantities discussed from here on depend on the time difference, t 



— t', but we nevertheless use the notation tt' for brevity. 
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with 

^^^^ "2^/ (2^) 

Note that P{t 0) = l/27r, i.e. / dujP{uj) = 1. 

In order to complete the argument, the phase fluctuations are easily related to the potential fluctuations, 
according to the relation (compare Eq. ( fTOl l) 

J«(0 = {^iixt)ip2{xO))^ = [D-lD-,],{Mxt)Mx'0))^h^,',t'=o , (21) 

and the potential fluctuation can be related to the interaction. This last step, however, requires some 
discussion. First, note that expanding the exponent in Eq. ^ up to second order, neglecting higher order 
terms, corresponds to the random phase approximation (RPA). Given this (Gaussian) approximation, the 
procedure for performing the phase average, outlined above, is permitted. As an important point, however, 
the RPA is known to be exact for calculating the density response of the Luttinger model, and hence the 
effective interaction, which concludes the argument: the potential fluctuations are given by 

{<^i(a;i)02(a;'t'))$ = '-V'^{xt,x't') (22) 

where V^{xt^ x't') is the screened (RPA) interaction. Thus — {i/2)D^^D^^V^ in an obvious short- 
hand notation. 



3 Effective impedance and bulk DoS 

The approach described in the preceding section is very similar to the theory of charge tunneling in ultra- 
small junctions, as reviewed, e.g., in Iil7i . which is also known as P{E) theory. (See also |,18J for related 
articles.) In this context, the quantity P{E) characterizes the influence of the environment on the tunneling 
between the two electrodes. For example, the forward tunneUng rate is found to be given by 



f{V) = [e^Rry^ / dEdE' f{E)[l - f{E' + eV)] P{E - E') 



(23) 



where V is the voltage, Rt the tunneling resistance, and f{E) the Fermi function. Note the detailed 
balance condition P[—E) = e~^^ P[E). Physically, P[E) describes the probability of exchanging the 
energy E with the environment ifTTl . 

On the other hand, in the present context, consider the tunneling rate from the interacting one-dimensional 
wire into a non-interacting lead, which clearly is given by 



T[V) = (e^i?^)"! / def{e)N{e)[l - f{e + eV)] 



(24) 



with N{e) the normalized DoS as given in Eq. iT% . Using this equation as well as the properties of P{uj) 
as described in the previous section and the detailed balance condition, it is straightforward to confirm that 
Eqs. ( l23T l and (l24l i coincide, provided P{E) from the tunneling theory is identified with P{u!)- The reason 
for the complete correspondence between these two quantities is apparent on physical grounds, since the 
quantity P{uj) characterizing the Luttinger liquid arises from the auxiliary potential fluctuations due to the 
fermion-fermion interaction. In the zero-temperature limit, one finds easily 



N{e) 



dujP{uj) (T = 0) 



(25) 



© 2003 WILEY- VCH Verlag GmbH & Co. KGaA. Weinheim 



5 



Exploiting the analogy further, we identify 



OJ 



Rk 



(26) 



where Rk = h/e^ = 27r?i/e^ is the Klitzing constant, and we may call Z{uj) effective impedance of the 
Luttinger liquid. For the simple example of an ohmic impedance with a high-frequency cut-off 



ReZ(w) 



1 



1 



Rk 

the result is ifTTl 



1 



T{2/g) u 



2/s 



{T ^0, <u <ujr) 



where 7 is Euler's constant. Thus the DoS is found to vanish at the Fermi energy, N{e) 
the exponent is given by 



9 



{f 



UJ 

-I < - 



dq 



TT J 27r {—iuj + + iqvp)'^ 



(27) 
(28) 

|e|2/9, where 
(29) 



where V^{qLLj) is the screened retarded interaction. In order to determine this quantity, consider the 2x2 
matrix structure of right- and left-moving particles, denoted by "+" and "— ", the bare interaction and the 
non-interacting response function: 



94: 92 

92 34 



Xo 



2'KVF 



—LJ — iO-irqVF 







qvF 



+uj-\-iO+qVF 



(30) 



as well as the RPA equation = (i + VoXo)^Vo. This matrix structure implies that we have to introduce 
a branch index for the fluctuating fields, which we suppressed up to now: ip,'b ^ (p±, <l>±. The density 
of states of the right-moving particles, for example, is determined from the (^+-(^+ correlation functions. 
The relevant quantity, V^_^, is straightforwardly determined and inserted into Eq. (|29] l: the q-integral can 
be done by contour integration, and we obtain the standard |fT9l result|3 



2/g^{K + K-^ -2)/2 



(31) 



4 DoS boundary exponent 

Near a boundary or an interface the quasiclassical approximation is insufficient, and the quasiclassical 
propagators pointing into or out of the boundary (or interface) have to be connected by boundary condi- 
tions ll2n |5| . The boundary condition is particularly simple for spinless fermions in one dimension at an 
impenetrable wall: obviously gu' {px) = gu' {—px), to ensure that there is no current through the wall. 

In the following we assume that the wall is located at a; = 0, and that the particles move in the half-space 
a: < 0. Instead of considering two branches of fermions in this half-space, we find it more convenient to 
mirror the left movers at the boundary and to consider only right movers in the full space, i.e. 

, ^ - j 9tt'{+,x) fovx<0 

[ gw{-,-x) forx>Q ^^2) 

+ ^^> \ $_(a;) fora; > ^ ^ 



^ See 1201 for a recent summary. Note that often the parameter K is denoted by g, which we avoid here since we prefer to use 
the latter symbol for the dimensionless conductance, see Eq. )27K in accordance with the P{E) theory 1 17 1. 
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This Green's function solves Eq. (|5]l both for x > and a; < 0. The bare interaction is now given in real 
space by 



Vo{x, x') = gA,5{x — x') + 52(5(2; + x') . 



(34) 



The Fourier transfrom of Vo is thus off-diagonal in momentum space, since the g2-term couples q with —q. 
Considering now the 2x2 matrix structure of particles with momentum q and —q, the bare interaction is 



(35) 



i.e. identical to what was given in Eq. ( [30l ) above with just a different meaning of the matrix index. Also 
the non-interacting response function and the RPA equation for the screended interaction are the same as 
before. Due to the absence of translation symmetry the effective impedance depends on the distance from 
the boundary, and is given by 









54 


92 1 




Va-q.-q , 




92 


9i J 



ReZ{x,uj) 
Rk 



LO 

2^ 
1 

4 



-Im 



dq Vj^J^) 



^2iqx _ 



yR 



(CO) 



27r {—iu! + iqvpY {—iuj + iqvF){—iuj — iqvp) 

[K + K-'' -2 + cos{2ux/v){K-^ -K)] . 



Hence the density of states at the boundary is found to vanish as 



N{e) 



\{1-K)/K 



(36) 
(37) 

(38) 



in agreement with the boundary exponent obtained in fTT\, and the impurity exponent given in fTS\. In 
fact this result has been considerably debated 120. .23. .24J : compare also Ii25il . as well as the detailed 
presentation by von Delft and Schoeller [2610 



5 Summary 

Quasiclassical theory is known to be a useful tool in the theoretical description of superconductors and 
superfluids, including non-equilibrium states and interfaces. For normal-conducting electrons the theory 
mainly serves for the microscopic foundation of a Boltzmann-like transport theory. However, when taking 
into account the Coulomb interaction in terms of a fluctuating Hubbard- Stratonovich field, the theory also 
captures important interaction effects, which are beyond the reach of the Boltzmann equation. In this 
article we illustrated this by considering the one-dimensional version of the theory, and demonstrated that 
it is possible to describe the non-Fermi liquid physics of Luttinger liquids. The mathematics involved is 
very similar to the P{E) theory of tunneling, which we made use of in this article, and also to the functional 
bosonization approach to Luttinger liquids [281 (which we did not exploit here). We concentrated on the 
spinless case and on thermal equilibrium; the spin is straightforwardly included in the theory, and clearly - 
since we use the Keldysh formalism - non-equilibrium situations can be covered as well. We are confident 
that in the future the quasiclassical theory will become a useful, complementary approach to transport in 
interacting one-dimensional systems. 
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